Quantum measurement in the charge representation 

J. Rammer, A. L. ShelankovQ and J. Wabnig 

Department of Physics, Umed University, SE-901 87 Umea 
(Dated: 21 December 2003) 

Counting statistics of charge transfers in a point contact interacting with an arbitrary quantum 
system is studied. The theory for the charge specific density matrix is developed, allowing the 
evaluation of the probability of the outcome of any joint measurement of the state of the quantum 
system and the transferred charge. Applying the method of charge projectors, the master equation 
for the charge specific density matrix is derived in the tunneling Hamiltonian model of the point 
contact. As an example, the theory is applied to a quantum measurement of a two-state system: 
The evolution of the charge specific density matrix in the presence of Nyquist or Schottky noise is 
studied and the conditions for the realization of a projective measurement are established. 
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I. INTRODUCTION 

In recent years, interest in quantum measurement has 
emerged in the context of solid state devices. This is in 
part forced by practical issues in connection with exper- 
iments involving nanostructures, and general questions 
of decoherencCj'^ and issues arising in nanomechanics, for 
example in probing the quantum electro-mechanical be- 
havior of a nano-resonator, and is also crucial for the 
read-out of a quantum computational process. Solid state 
nanodevices such as quantum dots are candidates for im- 
plementation of spin qubitSfSiMi^ and charge qubitSfS*i 
and superconducting nanodevices containing Josephson 
junctions have been studied in detailf^ and are being 
tested for their potential as charge qubitsiSiiSiii Nanome- 
chanics and monitoring qubits thus confront us with the 
practical details of a quantum measurementi^ In a quan- 
tum measurement one must, for the system purported to 
function as a detector, identify a collective variable which 
behaves classically. Such a variable is not a priori pro- 
vided by quantum theory, but in an electronic device the 
candidate is ultimately related to the charge flow. It is 
therefore of importance to have a quantum description 
of the charge transfer statistics, and demonstrate that a 
charge measurement can be a measurement of the state 
of a quantum object coupled to it. The purpose of the 
paper is to base such a description directly on the density 
matrix of the many-body system. 

The renewed interest in the quantum measurement 
process has in particular focussed on a two-level sys- 
tem coupled to a quantum point contact in the limit 
where the point contact can be modelled as a low trans- 
parency tunnel junction. The realization of the two-level 
system could be two coherently coupled quantum dots 
between which an excess electronic charge can tunnel.^ 
The two dots are electrostatically coupled to the tun- 
nel junction, or quantum point contact, thereby making 
the tunnelling amplitude depend on the state of the two- 
level system. This model has recently attracted much 
attentionpiiiiii^iiSiiLiSiia Gurvitai^ derived for the zero 
temperature case, a Markovian master equation for the 



density matrix of the two-level system keeping track of 
the charge transferred through the junction. Goan et. 
alm^ considered the Bloch-Redfield equation for the spin 
dynamics, i.e., the electron degrees of freedom of the tun- 
nel junction are traced out. To account for individual 
tunnelling events they employed the quantum jump ap- 
proach often used in quantum optics, and obtained from 
the stochastic master equation the spin evolution for spe- 
cific realizations of tunnelling events^ To account for 
specific realizations of tunnelling events, Korotkov sim- 
ilarly employed a stochastic treatment and showed by 
numerical simulation how an initially mixed spin state 
can evolve into a pure stateJS Rnskov and Korotkov 
considered the Markovian master equation for the den- 
sity matrix for the two-level system keeping track of the 
charge transferred through the junction at finite tem- 
peratures and calculated the noise spectrum>ii Recent 
achievements in counting statistics of charge transfers^ 
were used by Shnirman et. al.^ who derived a mas- 
ter equation for the spin density matrix, keeping track 
of the charges passing through the tunnel junction by 
a counting field as practised in counting statistics, and 
calculated the noise spectrum of a quantum point con- 
tact, with the same result as obtained by Bulaevskii et. 
alm^ Recently, a quantum oscillator coupled to the junc- 
tion has been studied, but by a method only valid at zero 
temperature.^^ The witnessed variety in techniques, calls 
for a standard method to attack these types of problems. 

We shall in the following develop a regular method 
for a joint description of the charge kinetics of a tunnel 
junction and a quantum system coupled to it. Essen- 
tial to the approach is that for an arbitrary many-body 
system we show how to treat the number of particles 
in a given spatial region or piece of material as a de- 
gree of freedom. This is achieved by employing suitably 
constructed charge projection operators introduced pre- 
viously in the context of counting statisticSf^S We shall 
refer to this reduced description, where at any moment 
in time the probability distribution for the number of 
particles or charges in a chosen region is specified, as the 
charge representation. The key construct of the theory 
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is the charge specific density matrix where the degrees 
of freedom of the environment are partially traced out. 
The charge specific density matrix allows one to evaluate 
the probability of the outcome of any joint measurement 
of the system and the charge state. Applying the charge 
projectors, allows us to use a standard kinetic approach 
to obtain the equation of motion for the density matrix 
in the charge representation. 

Tunnel junctions and quantum point contacts are ubiq- 
uitous in nanodevices such as for example those utilizing 
two-dimensional electron gases in semiconductor inver- 
sion layers combined with split gate technique. The de- 
veloped charge specific density matrix description of the 
dynamics of a quantum object we believe is the opti- 
mal description of such nanodevices since in electrical 
measurements any information beyond the charge distri- 
bution is irrelevant. This sets the stage for considering 
questions regarding quantum measurements, for exam- 
ple whether a measurement of the charge can provide a 
measurement of the state of a quantum object coupled 
to the junction. Indeed we shall demonstrate that the 
junction is able to perform a projective von Neumann 
measurement. The model allows us to study analytically 
a quantum measurement in the proper language of the 
density matrix. The intrinsic quantum bound on the 
measurement time necessary for the tunnel junction to 
operate as a measuring device is established. Quantum 
theory and measurement has a long and controversial his- 
tory, and the model provides another simple illustration 
that a quantum measurement can be described in stan- 
dard quantum mechanical terms, in fact in full detail for 
a realistic model of a nanodevice. The measuring scheme 
illustrates that there is no need to postulate the classi- 
callity of any variable, the charge state of the junction 
being directly accessible, nor to invoke "wave function 
collapse." The transmission of the electrical noise in the 
tunnel junction to the quantum object turns out to be 
sufhcient for the density matrix to decohere in the pointer 
basis. Amplification from the quantum to the classical 
level, the emergence of a projective measurement, can 
thus be followed in a realistic model of a nanodevice. 

The paper is organized as follows: in section II we 
construct the charge representation. In section III, the 
coupled object-junction Hamiltonian is introduced, and 
in section IV the equation of motion in the object and 
charge variables is obtained. In section V, we briefly dis- 
cuss an isolated tunnel junction in the charge represen- 
tation. The method we present is quite general, but for 
illustration we shall in this paper apply it to a low trans- 
parency tunnel junction coupled to a two-level system or 
for short a spin, and thereby obtain a master equation in 
terms of the relevant variables, charge and spin. In sec- 
tion VI, a short time measurement of the two-level system 
is considered and an analytic solution of the master equa- 
tion is obtained. Expressions for the characteristic times 
for decoherence and spin-charge separation are obtained. 
In section VII, the temporal progression of a quantum 
measurement is studied and the emergence of a projec- 



tive measurement is seen to be different depending on the 
relationship between voltage and temperature. In section 
VIII, we estimate the intrinsic quantum bounds on the 
measurement time. Finally, we summarize and conclude. 
Details of calculations are presented in appendices. 



II. CHARGE REPRESENTATION 

We start by showing that the number of particles in 
a given spatial region can be treated as a degree of free- 
dom. We construct the probability distribution for the 
number of particles in a specified region for an arbitrary 
many-body system. This is accomplished by using charge 
projection operators originally introduced in Ref. p^ . 

Consider an N-particle quantum system in a volume 
V and in state ^E". According to quantum mechanics, the 
probability, , that n particles are in a sub- volume A is 



.drN |^'(ri,...,rAr)|2 e^_^(ri, . . . , rjv) 

(2.1) 

where the indicator function 0^^(ri, . . . , r^v) is equal 
to 1 if exactly n of its N spatial arguments belong to 
the region A, and is zero otherwise. The volume A and 
its supplement volume B partitions the total volume into 
two non-overlapping parts, and the probability that there 
are N — n particles in region B equals the probability for 
having n particles in region A. 

The probability can be written as the expectation 
value 
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(2.2) 



of the Hermitian operator which acts in the position 
representation on the wave function according to 

V^-f{ri, . . . ,rAr) e,f_jv(n, . . . ,rAr)^'(ri, . . . , r^). 

(2.3) 

From their definition, one observes that the introduced 
operators, V^, have the properties 
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(2.4) 



and therefore constitute a complete set of Hermitian pro- 
jectors. The projected state \^n) = T^nl^) is the com- 
ponent of 1 5*) for which exactly n of the particles are in 
region A. 

The many-particle operator can be expressed in a 
simple way through single particle operators. We intro- 
duce the gauge transformation operators 



Ux = exp 



N 



k=l 



(2.5) 



where A is a real parameter, and 6^{rk) equals 1 if 
belongs to the region A and is zero otherwise. By virtue 
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of the identity, 

/ = e;^,^ (2.6) 



the projection operator can be presented in any repre- 
sentation as 

27r 

^" = / ' (2.7) 



thus presenting the many-particle operator in terms 
of the product of the single particle gauge transforma- 
tions in U^. This result can also be obtained by observ- 
ing that the projected state, [vPn), is an eigenstate of U^, 
U^l'^n) — e''^"|5'„). If the region A, where the particles 
are counted, occupies the whole volume V, the charge 
projector coincides with the one used by P.W. Ander- 
son in the theory of superconductivity for projecting the 
BCS-wave function onto the state with a fixed particle 
number. 



A. The charge specific density matrix 

We shall quite generally consider a system consisting 
of two parts: a quantum object to be measured and the 
measuring device, which is taken to be a tunnel junction 
connecting two electron reservoirs. 

For the sake of the derivation, we assume that the 
object together with the electron reservoirs constitute a 
closed quantum system. The system is described by its 
full density matrix p{£^, Rei]^' , R'^i), where ^ is the co- 
ordinate of the object and Rgi = ri,...rff comprises 
the electron coordinates, N being the total number of 
electrons relevant for the functioning of the measuring 
device. With the understanding that the only accessible 
information about the electrons is their charge distribu- 
tion, we introduce the charge specific density matrix of 
the measured system 

Pn = Trel {Vn p) (2.8) 

where Vn is the projection operator introduced in the 
previous Section, and the trace is with respect to the elec- 
tron degrees of freedom Rei . The charge specific density 
matrix p„ = PniS.,0 allows one to deduce the probabili- 
ties of any joint measurement performed on the quantum 
object simultaneously with a charge measurement. 

If the charge specific density matrix is traced over the 
remaining degrees of freedom, the probability Pn that 
there are n charges in a specified region (reference to 
which is suppressed in the following) is the expectation 
value of the charge projector, or expressed in terms of 
the charge specific density matrix 



where the trace is with respect to the degrees of free- 
dom of the quantum object, i.e., ^. It follows from 
Eq. (|2.4|l that the probability distribution is normal- 
ized, X^nP" ~ I- Apart from proper normalization, /5„ is 
the density matrix for the quantum object after a charge 
measurement. 

The charge projection operators define what we shall 
refer to as the charge representation. The degree of free- 
dom in the charge representation is the variable n: the 
number of charges in a specified spatial region or equiva- 
lently in a piece of material. All other information about 
the charges is traced out. The charge representation is 
a strongly reduced representation, the description of the 
environment is reduced to one variable, but if interest 
is in the currents in a system, information beyond the 
charge representation is irrelevant. 

We now consider the charge specific dynamics, i.e., the 
equation of motion of the quantum object for specified 
charge. The dynamics of the combined system is gov- 
erned by the equation for the density matrix 

ip=[H,p] (2.10) 

and the equation of motion for the charge specific density 
matrix is thus 

i'p„=Tr,i {Vn[H,p]). (2.11) 

The charge specific density matrix describes the quan- 
tum dynamics of the coupled object conditioned on the 
charge variable being specified. Since the charge specific 
description is a reduced description, generally a hierar- 
chy of equations is generated. In the following we shall 
study a situation where /5„ can be obtained explicitly. 

III. MODEL HAMILTONIAN 

In this section we discuss the Hamiltonian for an arbi- 
trary quantum system coupled to a tunnel junction 

H = Hs + Hi+Hr + HT. (3.1) 

Here, Hg is the Hamiltonian for the isolated quantum ob- 
ject, the system to be measured, and a hat marks terms 
which are operators with respect to the degrees of free- 
dom of the quantum object. The Hamiltonians Hi^r spec- 
ify the isolated left and right electrodes of the junction 

Hi=^eiclci , Hr = ^erclcr (3.2) 

I r 

where Z, r labels the quantum numbers of the single par- 
ticle energy eigenstates in the left and right electrodes, 
respectively, with corresponding energies ej,^ and annihi- 
lation and creation operators. The operator Ht describes 
tunnelling. 



Pn = Tr^ (pn) 



(2.9) 



(3.3) 
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where 



IrCj C. 
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(3.4) 



inl- 



and are operators acting on the coordinates, ^, of the 
quantum object. The Hermitian property of the Hamil- 
tonian requires that Tir — T^i- Compared to an isolated 
tunnel junction, the additional feature of the model is 
that the tunnelling amplitude, and thereby the conduc- 
tance of the tunnel junction, depends on the state of the 
measured system. 

In the next section we shall study the dynamics of the 
system when the probability distribution for the number 
of charge transfers through the junction is specified at 
all times. This is achieved by using the charge projectors 
introduced in section^ For a tunnel junction, the spa- 
tial region of interest for counting charges is either of the 
two electrodes, say we choose the left one. Of importance 
are in view of Eq. (|2.11|) the commutation relations of the 
charge projectors and the Hamiltonian. The terms in the 
Hamiltonian commute with the charge projectors except 
for the tunnelling term. The discrete charge dynamics of 
the tunnel junction is specified by the charge projection 
operators according to 



VnClCr ^ clcrVn^l , :P„4q = 4qP„+i. (3.5) 

In terms of the tunnelling operators, Eq. (|3.4() . the iden- 
tities read 



(3.6) 



These identities are used repeatedly in the derivation of 
the equation of motion in the charge representation. 



IV. CHARGE SPECIFIC DYNAMICS 

In this section we shall obtain the equation of motion 
for the charge specific density matrix, p„, to lowest order 
in the tunnelling. 

Taking advantage of the relations in Eq. (|3.5|) . the 
equation of motion for the charge specific density ma- 
trix, Eq. H2.11|l . can be written in the form 

Kit) + i[Hs,p„{t)] 

= E {TlrA[:\t) + frlBi-\t)) + h.C. (4.1) 
Ir 

where here and in the following h.c. represents the her- 
mitian conjugate term, and the time dependent system 
operators A and B are given by 



(4.2) 
(4.3) 



One sees from Eq. 1)4. that the time evolution of the 
charge diagonal component of the density matrix, pn, is 
determined (in addition to the internal dynamics of the 
system) by the charge off-diagonal components Vn±ip'P 
that control the matrices A and B. The latter two are 
small, being generated by the rare tunnelling events, and 
can be expressed in terms of the diagonal elements using 
perturbation theory. To obtain, e.g. A, one uses the 
equation of motion, Eq. H2.10|l . to obtain 

A,^ it) - ^c.,,A[:\t)+^[Hs,A[:\t)]^ 

- Tt,i (clcr[HT,p{t)]rn) (4.4) 

where ioir = ei — £r- To lowest order in tunnelling, one 
evaluates the source term on the r.h.s. retaining only 
the charge diagonal VmP'Pm components (m = n and 
m = n — 1) of the full density matrix p; the off-diagonal 
components, Vn-2pVn and Pn-ipVn+i, give higher or- 
der corrections which are neglected (for details we refer 
to appendix . The diagonal terms can be expressed 
in terms of the single particle distribution functions for 
the electrodes, fi^r, which, since the electrodes act as 
particle reservoirs, can be taken independent of the total 
charge number of the particles in the left electrode, n. As 
shown in appendix to lowest order in the tunnelling 
the inhomogeneous term on the right in Eq. (|4.4(l can be 
expressed in terms of the charge specific density matrix 
and an explicit solution for A (and similarly for B) can 
be obtained. The resulting equation is a non-Markovian 
master equation for the charge specific density matrix for 
the system, Eq. (|A15p . In the limit where the part of the 
tunnelling matrix element in Eq. H6.1|l which depends on 
the coupling to the system is small, the temporal non- 
locality of the kernels can be neglected and the Marko- 
vian equation is obtained. The master equation for the 
charge specific density matrix /5„ has the form: 

Kit) = ^[iHs + M),Pnit)] 

+ C{pnit)}+V{p:{t)}+J{p'jt)} (4.5) 
where and p'^ denote "discrete derivatives" , 



Pn — 2 iPn+1 ~ Pn~l) , 
Pn = Pn+l + Pn-1 - '^Pn- 



(4.6) 
(4.7) 



The Hamiltonian Hs describes the dynamics of the iso- 
lated system, and the rest of the terms on the r.h.s. of 
Eq. H4.5|l describe the effect of the tunnelling events; in 
the formulae below, the subscript +(— ) marks the contri- 
bution of the processes when an electron tunnels from the 
left (right) electrode to the right (left) one. The "mag- 
netization" term, M = M, + M_, 
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T, 



Ir 



Ir 



Tl 



- h.c.j (4.8) 

- h.c) (4.9) 
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renormalizes the system Hamiltonian Hs- The operator C = + where C± {p} maps p according to 
I 



Tl 



+ h.c. 



(4.10) 



and 



Ir 



ft 



h.c 



(4.11) 



The diffusion, and drift, operators are 
j{p]=V+{p]-V-{p], 



(4.12) 



where 



^+ = E /'(I - f^/rlpT'i- + h.c. (4.13) 



and 



=^/,(l-/j)Ti, 



Tl 



(4.14) 



ir 



We have introduced the notation (but suppressed the 
time dependence in the above formulas) 



ft 



it) 



dr e 



icJi-r-T+i j dt'V{t') 



(4.15) 

and assumed a voltage U applied to the junction, V{t) = 
eU{t), e being the electron charge. 

The master equation, Eq. (|4.5|l . requires an initial con- 
dition. We assume that the left and right electrodes of 
the junction are disconnected at times preceding the ini- 
tial moment i = 0, so that the electrodes are in definite 
charge states, say with n^^^ electrons in the left electrode. 
Using the convention that n is counted from n^^\ the 
variable n becomes the number of electrons transferred 
from the right to the left electrode. The initial condition 
reads 



Pnit = 0) = Sn,0 P 



(0) 



(4.16) 



where p^*^^ is the initial density matrix of the measured 
system. 

The kernel V describes diffusion in charge space, and 
the average charge current from the left to the right elec- 
trode, — e^nTr /5„(i), is given by the drift term J 



I{t) = eTr J {m} 



(4.17) 



where p is the (unconditional) density matrix 



(4.18) 



We stress that the form of the derived master equa- 
tion for the charge specific density matrix is valid for 
any quantum object coupled to the junction. In the 
present paper, we shall consider the electrons cou- 
pled to a two-level system or for short a spin, a sys- 
tem considered in numerous papers as discussed in the 
introduction,iiiiiiSiiSiiLiS*ia and /5„ will be a 2 x 2 ma- 
trix. A remaining task is thus to unravel the complex 
spin structure of the master equation presently hidden 
in the bracket operation. However, first we consider the 
isolated timnel junction. 



V. ISOLATED TUNNEL JUNCTION 

To familiarize with the charge representation we pause 
to consider the case where the system is decoupled from 
the junction. In this case, the only degree of freedom of 
interest is the charge n, and the equation for the proba- 
bility, pn, for the transfer of n charges 



Pn 



DPn + IPn 



(5.1) 



is obtained by taking trace of the master equation 
Eq. H4.5|l . Here, 



I^2TTj2\Tlr\Hfl-fr)S{ei+V- 



and 



D ^7rJ2\Tlr\^ifl +fr- 2flfr)5{ei + V - 



(5.2) 



(5.3) 



For simplicity, we assume the bias V is time independent. 

The physical meaning of the parameters / and D in 
Eq. (|5.1|l can be deduced from the structure of the equa- 
tion. Multiplying Eq. H5.1|l by n and performing summa- 
tion with respect to n, one readily sees that / equals the 
dc current through the junction, i.e., the average rate 
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of charge transfers from the left to the right electrode, 
I ~ —-^Yl ""^Pnit)- The value of the dc current / — I{V) 

n 

is given by Eq. H5.2|l . For the time evolution of the vari- 
ance of the charge distribution, 

An{t) = {n^t)) - {nit)y , (5.4) 

where {n'^{t)) — '^n^Pnit), one obtains from Eq. (jS.lll 



dt 



2D 



(5.5) 



thus D has the meaning of the charge diffusion coefficient, 
characterizing the randomness of the charge transfers. 

To connect with the standard noise discussion in terms 
of current fluctuations, we express the transferred charge, 
n{t) as the time integral of the current i(t), n{t) = 
J* dt'i{t'), and the variance takes the form 



t' 



An(t) = 2 j dt' j drS^T) 



(5.6) 





where S{t) = S{—t) is the current-current correlator, 

S{t) ^ {6i{t)5i{t + t)) , 5i{t)^i{t)- {i). (5.7) 

At times much larger than the current-current correlation 
time, one obtains from Eq. H5.6|l 



dt 



OQ 

2 J drSir) 



(5.8) 



Recalling Eq. (|5.5|) . one establishes the relation, 

4.D = S^^o , (5.9) 

between the charge diffusion coefficient D 
and the power spectrum of the current noise, 
Sui = 2 drS^r) cos{lut), at zero frequency. 

In quasi-equilibrium, when fi^r are Fermi functions 
corresponding to the temperature T, one readily ob- 
tains from Eqs. H5.2() and H5.3|l that the charge diffu- 
sion constant D is related to the dc current / as 2D — 
Jcoth V/2T. In view of Eq. (|5.9|l . this corresponds to the 
well-known result 



S-^^o = 2/(F) coth ^ 



(5.10) 



which expresses the fluctuation-dissipation theorem (and 
holds for arbitrary voltages for the case of a tunnel junc- 
tion )^ 

Given the initial charge distribution, Eq. (j5.1|l can 
be solved introducing the Fourier transform x(A;t) = 
J2nPn{t)^^'^^ ■ For the initial condition p„(t = 0) = Sn.o, 
the solution reads 

x(A;t) = exp[2Di(cos A - 1) - ilt sin X] . (5.11) 



Not surprisingly, this expression reproduces the long time 
(Markovian) limit^^ of the generating function of the full 
counting statistics for a tunneling junction first derived 
in Rcf. 26. Inverting the Fourier transform (see appendix 
IB II for details), the probability of n-charge transfers in 
the time span t reads 



/ /n) \ 

Pn{t) = /„ ^- — exp [vn - (n) cothu] (5. 
\ smh V J 



12) 



where /„ is the modified Bessel function, (n) = I{V)t is 
the average charge transfer for the bias V, and v = V/2T. 
The moments of the distribution, (n'^) and (n— (n))'^, can 
be evaluated with the help of the generating functions 
Eq. IjBQp in appendix IB II 



VI. SHORT TIME MEASUREMENT 

The charge specific density matrix can be used to de- 
scribe a quantum measurement in a two-fold way: (i) 
measuring the charge transferred through the junction 
provides information about the state of the system cou- 
pled to the junction, or (ii) a measurement of the coupled 
system can reveal information about the charge state of 
the junction. In this paper we shall consider the case 
where the system coupled to the junction is a two-level 
system or for short a spin, whose internal dynamics is 
specified by Hg = ila^ ■ We are interested in the read-out 
of the state of the spin, and we analyze how measuring 
the charge can be a projective measurement of the spin. 

For definiteness, we chose to measure the z-projection 
of the spin, and for the measurement in this basis the 
coupling to the spin is taken to be 



Tlr 



Vlr + CTzWlr 



(6.1) 



so that the spin up and down states correspond to ex- 
treme values in tunnelling strength. It is convenient to 
introduce the following quantities: 



r Gv ' 






1 Gw 
1 




Wlr 
Uur ' 


[ G2 J 


Ir 


. U2lr . 



dfjei) 
dsi 



S{si - er) 



(6.2) 

where Vir = |wjr-P, Wir = |wjrP, Uiir = ^Vi^wir, and 
U2ir = '^Vi^wir- For the up and down spin orientation, 
the tunneling conductance (in the units /K) is given by 
Gv + Gw ± 26*1 , respectively. 

Following von Neumann,^'' we aim at an effectively 
instantaneous measurement of the spin quantum state, 
which means that the time of the measurement must be 
short compared to the intrinsic spin precession time. To 
describe the measurement process, we are therefore in- 
terested in solving the master equation for times much 
shorter than the inverse of the Rabi frequency fi. In that 
case, Hs in Eq. H4.5|l and Eq. (|4.15|) can be set to zero, 
and the spin structure of the master equation, Eq. (|4.5|l . 
reduces considerably. 
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In this limit, the renormahzation term Ai = M-+ + 
M_, defined by Eqs. H4.8|) . and (|4.9|) is = M^dz where 

fi - fr 



Wlr 



Ir 



El + V — Er 



(6.3) 



The physical mechanism behind the renormahzation of 
the spin "Zeeman" energy given by is the sensitivity 
of the total electron energy to the spin orientation. In- 
deed, in second order perturbation theory with respect to 
tunneling, the single electron levels acquire a shift con- 
trolled by the square modulus of the tunneling matrix el- 
ement which is either \vir + P, when the spin is up, or 
\vir — wir\^ , when the spin is down. The interference term 
2^Vi^wir gives rise to the spin dependence of the energy 
and, therefore, in Eq. (|5.3|l . The sign of the combina- 
tion ^Vi^wir is not fixed by any physical requirement. In 
a mesoscopic system, we expect it to be a random func- 
tion of the quantum numbers I and r, fluctuating both 
in absolute value and sign. In view of this, the contri- 
butions to Mz from the states in the wide energy range 
of the order of the Fermi energy, have the tendency for 
cancellation. As an estimate, ~ GiSE where SE is 
a narrow energy interval within which ^Vi^wir are cor- 
related. The renormahzation may therefore be a small 
correction, which we neglect below. 

The master equation, Eq. (|4.5|) . becomes 

G2V V 

Pn = — [^z,Pn] + 2GwV COth — {a^Pn^z ~ Pn) 

I 21 

+ ^Vcoth ^ (GvPn + Gi[/5::,CT,]+ + Gw^zPn^z) 
+ V {GvPn + Gl[Pn7^z]+ + Gw^zPn^z) 

+ iG2V[p:,az]+iG2V COth ^[p[„az], (6.4) 

where Gv,w,i.'^ defined in Eq. (|6.2|1 . 

At sliort times, the components of the density matrix, 



(6.5) 



are decoupled, and solving Eq. (|6.4|l amounts to solving 
the following three equations: 



iin = lvcoth-^{Gv + Gw + 2Gi)u';^ 



2 2T 
V{Gv + Gw + 2Gi)u'^ 



1 V 
-l^coth — (Gv 

2 2T^ 



+ V{Gv + G 



Gw — 2Gi)d'^ 
w - '2Gi)d'n , 

V 



(6.6) 



(6.7) 



-2iG2Van - 2GwVcoth— q;„ 
i (vcoth^iGv-Gw)-2iG2V]a: 



V{Gv - Gw) - 2zG2l^coth ^ ) a'„. (6. 



Solving the master equation, Eq. I|6.4|l . is thus 
achieved once the equation with the structure 



(6.9) 



is solved. According to appendix^ the solution is spec- 
ified in terms of the modified Bessel function 

Xn{t)^xJ-^\ e-2^==*/„(2iv'-D:.+-D.-) (6.10) 

where Dx± — ± i J^,. All moments of the stochastic 
process a;„ (t) can be expressed in terms of the parameters 
Dx and which for the original problem becomes in 
terms of voltage, temperature and time. 



A. Characteristic times 

Let us now analyze the time dependence of the charge 
specific density matrix. It is found solving the master 
equation with the initial condition Eq. H4.16|> . where p'*'' 
is the initial spin state (to be measured). 



al do 



(6.11) 



First, we consider the unconditional spin density ma- 
trix p, Eq. H4.18|l . which gives information about the spin 
state irrespective of the outcome of the charge measure- 
ment (i.e., discarding the results of a charge measure- 
ment). The equation for p, 

G2V V 
P = — — i^z, p] + 2GwV coth — i&zpaz - p) , (6.12) 
I 21 

is obtained by summation with respect to n in Eq. (|6.4|l . 
Its solution reads. 



m 



Uq 



2iG2Vt -t/ra 



^^^-2^G2Vt^-t/r, 



where the parameter Td, 



tT^ = 2Gwycoth 



V_ 

2T 



(6.13) 



(6.14) 



is the decoherence time, which gives the decay rate of the 
charge unconditional off-diagonal elements of the density 
matrix J^n'^ri- 

As time passes, the probability distributions for being 
in the spin up or down states, Un(t) and dn{t), will ac- 
cording to Eq. I|6.10|) drift and broaden in charge space. 
To investigate the different drift in the two distributions 
we consider their overlap in charge space, '^Un(t)dn(t) 



and we get according to Eq. (|B4|I 
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where t = 2Ft coth ^^(Gy + Gw)^ - 46^ tanh^ ^. 

At large times, t ^ 1, the overlap decays exponentially, 
e~*/^== , with the characteristic rate 



= SV^tanh — 



2T Gl 



G 



w 




AGf tanh^ 



2T 



(Gv + GwY 



(6.16) 



After the time span Tz, the probability distributions for 
spin up and down have separated in charge spacei^ 

A similar calculation for the off-diagonal elements of 
the charge specific spin density matrix gives 



(6.17) 

where t* = 2tV^coth ^ |Gy — Gw + 2iG2|, and at large 
times the off-diagonal elements have decayed exponen- 
tially at the characteristic decoherence rate 



Gw — \Gv ~ Gw -\- 2iG2\) ■ 
(6.18) 

After a time span tj^ the coupling to the current has thus 
reduced the spin state to a mixture, the charge specific 
density matrix being diagonal in the measurement basis. 
We note that there is no simple relation between t± and 
the decoherence time Td in Eq. (|6.14|l . Generally, a„ is 
more robust relative to decoherence than a„ . 

Although the characteristic times in Eqs. H6.16|l . and 
(|6.18|) depend on different combinations of parameters of 
the model, their ratio seems to be rather universal and 
we expect quite generally that 



V 



— = c tanh — 
2T 



(6.19) 



where c is a constant of order unity»2S 

We have identified the characteristic time scales in the 
problem, and found their dependences on voltage and 
temperature. The three time scales describe the coherent 
dynamics of the spin under different experimental condi- 
tions. When the charge state of the junction is left unob- 
served, the dynamics of the spin is given by the uncondi- 
tional spin density matrix, Eq. 16.13|l . which decoheres on 
the time scale, r^, given in Eq. H6.14|l . When the spin evo- 
lution is conditioned on the charge state, the dynamics 
of the spin is given by the charge specific density matrix. 
In the shot noise regime, F 3> T, spin-charge separation 
and the decoherence of the charge specific density ma- 
trix happens on the same time scale, rj^ ~ Tz, and the 
tunnel junction provides a projective measurement of the 
spin after a time span of this order. For small voltages, 
y T, when <SC t^, the tunnel junction provides a 
projective measurement of the spin only after the long 
time span r^. 

To study the time evolution of the charge specific spin 
density matrix in detail we turn to analyze the obtained 
analytical results. 



VII. "WATCHING" A QUANTUM 
MEASUREMENT 



We shall now consider the time evolution of a measure- 
ment of the state of the spin performed by the tunnel 
junction. The charge specific density matrix is presented 
in the Pauli basis 



Pn{t) 



Pn{t) 



1 



(7.1) 



where Pn{t) is the probability that n electrons have 
passed through the junction at time t and the polariza- 
tion or Bloch vector s„(t) indicates a point within the 
unit Bloch sphere. 

The coupling between the spin and the tunnel junction 
is chosen to have the spin structure given in Eq. (|6.1|1 . 
and the parameters vir and wir are in the following taken 
to be real constants. In that case G2 vanishes; the ratio 
of the conductivities is taken to be Gw = O.OlGy. 

Initially the junction is in the definite charge state, 
Pn(0) = 6n,o, as depicted in figuren(d). For definiteness, 
we choose the spin initially in the pure state correspond- 
ing the Bloch vector So(0) = (—1,0, 0). The further evo- 
lution of Bloch vectors takes place only in the xz-plane, 
and the charge specific spin density matrix can there- 
fore be represented by a line rising vertically from the 
point indicating the location of the Bloch vector s„(i). 
Its height (ending in a black dot for visual clarity) repre- 
sents the probability for the transfer of n charges, Pn{t). 
The initial state is thus represented as depicted in figure 
n(a). At a later time t, many vertical lines are present, 
representing the probabilities for various charge transfers 
and their corresponding spin states as depicted in figure^ 
(b). In fact, many Bloch vectors s„(t) will appear close to 
each other, and placing a square grid over the xz-plane 
cross section several Bloch vectors, differing from each 
other in charge number by units of one, can lie in the 
same square. For visual clarity, however, only one ver- 
tical line in the figures will be associated with a square, 
and its height is the sum of the individual probabilities 
Pn belonging to the Bloch vectors in the square. 

We now turn to study the time evolution of the charge 
specific spin density matrix, and start by considering the 
shot noise regime where the junction is biased by a high 
voltage, y > r. 



A. Schottky limit 

The evolution of the Bloch spin density for the case 
where the voltage is larger than the temperature, V = 
20T, is presented in figure ^ (b) and (c). At time 
t — Tz/2 ~ t_l/2, the Bloch spin density is spread along 
the unit circle as depicted in figure^ (b) (the correspond- 
ing charge numbers increasing in the clockwise direction). 
The visualization shows that during the evolution, the 
charge specific density matrix stays pure. This is charac- 
teristic of the shot noise regime, where the charge trans- 
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fer statistics is Poissonian, though not universal since to 
some extent dependent on the choice of constant tun- 
nelhng matrix elements. At this time, t = Tz/2, where 
on the average roughly ten electrons have tunnelled, the 
charge probability distribution, p„ = m„ + c?„, depicted 
in figure is built equally of the charge conditioned 

probabilities for the spin to be in state up or down. At a 
larger time, t = Sr^, the Bloch spin density distribution, 
depicted in figure ^ (c) , has split and is located near ei- 
ther of the Bloch vectors S| = (0, 0, 1) or S| = (0, 0, — 1), 
corresponding to the spin up and down states, respec- 
tively. At this time, where on the average roughly one 
hundred electrons have tunnelled, the charge probabil- 
ity distribution, p„ = m„ + dn, is spin separated into a 
charge probability peak built solely by the probabilities 
for the spin to be in state up and a charge probability 
peak built by the probabilities for the spin to be in state 
down, and has the shape depicted in figure ^(f)- The 
total probability in the two peaks both equal one half, 
the probabilities initially for the spin to be in state up 
or down. The probability distributions for charge and 
spin have thus through interaction come in one-to-one 
correspondence, and measuring the charge state of the 
junction at times larger than is thus a measurement 
of the state of the spin prior to interaction with the tun- 
nel junction. The relative frequency with which a charge 
state in either of the two peaks is realized is equal to the 
probability for the corresponding spin state at the start 
of the measurement. The junction thus functions as a 
projective measuring device of the spin. To be effective, 
we see from figure ^ (f), that on the average hundred 
electrons must have tunnelled. 



B. Nyquist limit 

Next we study the evolution of the charge specific spin 
density matrix for the case where the voltage is smaller 
than the temperature, V — 0.2T. The initial spin and 
charge states are chosen as in the high voltage case de- 
picted in figure ^ (a) and (d), so initially the Bloch 
spin density is concentrated near the initial Bloch vec- 
tor, s — (—1,0,0). However at larger times, the evolu- 
tion of the Bloch spin density is quite different from the 
high voltage case proceeding in two steps as noticed from 
figure 121 (a). At time t±, the x-component of the Bloch 
vectors, s„(t_l), have started to decay, and the charge 
specific spin density matrix is no longer pure. At this 
time on the average only roughly five electrons have tun- 
nelled according to figure|21(d). At the much larger time 
t = Tz/2, Tz ~ IOOtj^, it is seen from figure|2](b), that the 
Bloch vectors are concentrated on the z-axis, the charge 
specific spin density matrix has evolved into a mixture. 
At this time, where on the average roughly one hundred 
electrons have tunnelled, the charge specific spin up and 
down distributions, m„ and dn, have started to charge 
separate, but not enough to significantly distort the to- 
tal charge distribution as evident from figure [21(e). At 



the later time t = 5tz, the Bloch spin density is again sep- 
arated and located at the sites corresponding to the spin 
up and down states as depicted in figure El (c) , and the 
charge probability distribution is again spin separated, 
and the total probability in the two peaks equals one 
half, the probabilities initially for the spin to be in state 
up or down. A projective spin measurement has again 
been performed by the tunnel junction. To be effective, 
we see from figure (21 (f), that on the average roughly one 
thousand electrons must have tunnelled. 

The evolution of the initial state is seen to take place 
in two steps, governed by the time scales t± and Tz. 
First the pure initial state, s = (—1,0,0), decays into 
a mixture where the Bloch vectors have no x-component, 
s = (0, 0, s^), (no off-diagonal elements in the charge spe- 
cific density matrix) . The mixture is then purified in the 
further time evolution, before the density matrix finally 
approaches the pure states |t) and ||) on the time scale 

Tz- 

VIII. MEASUREMENT TIME 

In the measuring scheme under consideration, infor- 
mation about the spin state is obtained from a charge 
measurement. In the long time limit, the charge prob- 
ability distribution acquires two distinct peaks (as seen 
in Figures n(f) and 121(f)), and the outcome of a charge 
measurement is almost always in the vicinity of or 
Ql^ Qld = being the positions of the peaks. Ide- 

ally, one reads off the z-projection of the spin, t or |, from 
the charge outcome Qt or Q^. At any finite time, this 
scheme may produce errors. The intrinsic error mech- 
anism is electrical noise due to which the peaks have 
a finite width and, therefore, partially overlap. Addi- 
tionally, the spin state is not pure at finite times, spin 
purification occurring only asymptotically. Also, errors 
may occur due to a finite resolution AQjj of the charge 
detector. The time span, tm, needed to make a reliable 
measurement, is defined by the condition that the infor- 
mation gained from the charge measurement suffices to 
identify the z-projection of the spin after the measure- 
ment. First, we shall consider the case of an ideal charge 
detector, before taking into account the finite resolution. 

a. An ideal charge detector, generates certain dis- 
crete values, n, as the outcome of the charge measure- 
ment with the probability p„ ~ Tr /5„, and leaves the 
spin in the state described by the density matrix pro- 
portional to pn- With pn of the form in Eq. H6.5I) . the 
probability p„ equals Pn = w„ + dn , and 

tin ~r i 

are the probabilities for the spin to be in the up or down 
states for a given observed n. 

Asymptotically, the distributions u„(i) and dn{t) do 
not overlap in accordance with Eq. Ht).15|l . Therefore, 
only one of the diagonal elements, w„ or c?„, remains finite 
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(a) t = 



(b) t = 



(c) t = 5tz 
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(d) t = 



(e) t = T^/2 



(f) t = 5r,/ 



FIG. 1: Spin-charge evolution in the high voltage case, V — 20T. (a) Spin-charge representation of the initial state. The 
arrow indicates the evolution of the unconditional spin density matrix described by the Bloch vector shrinking along the a;- axis 
decohering to the mixture state s = (0, 0, 0) on time scale Td. (b) As time passes the Bloch spin density spreads along the unit 
circle, (c) At time t = 5tz the Bloch spin density is located at the pure states ||) and |J.). Simultaneously, the initial charge 
distribution, depicted in figure (d), will due to tunnelling events start to move in charge space with a peak determined by the 
average number of charges transferred, (e) At time t = Tz/2, the charge distribution, Pn = u„ + dn, is built equally of the 
charge conditioned probabilities for the spin to be in state up or down, (f) At time t — dr^, the charge distribution is spin 
separated into a peak built solely by the charge conditioned probabilities for the spin to be in state up and a charge probability 
peak build by the probabilities for the spin to be in state down. A projective spin measurement has been performed by the 
tunnel junction after an average of roughly hundred electrons have tunnelled. 



and the spin is left in either the up or down state. At 
long times, the measurement is perfectly projective. 

To quantify the efficiency of the measurement at finite 
times, we take an approach well known in the context of 
information theory (see e.g. ref. l30f) . Considering the 
outcomes of the joint measurement of the charge n and 
spin a =t, i, as random variables, the probability for the 
outcome (n, "f) (or (n, J.)) is given by u„ (or dn)- The 
joint entropy, H{a,n)^ equals 



and the entropy H(n) of the charge variable n is, 

H{n) = -^Pn log2 Pn ■ (8.3) 

n 

The conditional entropy, H{a\n) — H{a,n) — H{n), 
is a measure of the information gained about the spin 
from a charge measurement. Combining Eq. (|8.2() and 
Eq. I|8.3|l . we obtain 

^^((^1") = -^Pn(pIlog2Pl+pilog2Pi) ■ (8.4) 



H{a,n) = (w„log2 



.log2d„), (8.2) 



This quantity changes from the initial value, 
— (wq log2 uq -l- do log2 ) 5 to zero at large times 
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(a) i = rx (b) t = ~ 50t^ (c) t = 5t^ 
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n n n 



(d) t = Tj_ (e) t = ~ 50t_i (f) t = 5t^ 

FIG. 2; (a) Spin-charge evolution in the low voltage case, V = 0.2T. (a) At time t = t±, the x-component of the Bloch vectors 
have started to decay, and the charge specific spin density matrix is no longer pure, (b) At the much larger time t = 

~ lOOrx, the Bloch vectors are concentrated on the z-axis, and the charge specific spin density matrix has evolved into a 
mixture, (c) At time t — 5tz, the Bloch spin density is located at the pure states If) and ||). (d) At time t = t±, only an 
average of hundred charges has tunnelled, and the charge distribution, p„ — Un + dn, is built equally of the charge conditioned 
probabilities for the spin to be in state up or down, (e) At time t = rz/2 ~ 50t±, the charge specific spin up and down 
distributions, u„ and d„, have started to charge separate, but not enough to significantly distort the total charge distribution, 
(f) At time t — 5tz, the charge distribution is spin separated into a peak built solely by the charge conditioned probabilities 
for the spin to be in state up and a charge probability peak built by the probabilities for the spin to be in state down. A 
projective spin measurement has been performed by the tunnel junction after an average of roughly one thousand electrons 
have tunnelled. 



where pjj"'- is either zero or 1. At intermediate times 
it is a measure of the error in the predicted spin state 
after the charge measurement. The time dependence of 
the conditional entropy for the initial state of the spin 
in the a;-direction is shown in Fig. |31 The curve shows 
exponential-like decrease of the conditional entropy from 
the initial value, one bit, to the value zero. 

For the case of an ideal charge detector, we define the 
measurement time tm as the time when the spin entropy 
decreases below a chosen threshold determined by the 
required fidelity of the measurement. In Fig. 01 we show 
the measurement time for various parameters with the 
entropy threshold (arbitrarily) chosen as 0.01. The main 
conclusion is that the intrinsic measurement time is of 
order as given in Eq. (|6.16|) . tm Tz- 

The most favorable condition for the accuracy of the 
measurement is the shot noise regime where the charge 



noise is relatively small. As a rough estimate of the mea- 
suring time, t,n T^, in the shot noise regime we obtain 
from Eq. 



7- 



eU 



8G? 



Gv + G 



w 



(8.5) 



where U is the voltage applied to the junction and 7 is a 
dimensionless constant. For the tunnelling Hamiltonian 
approach and the Markovian approximation to be justi- 
fied, the constant 7 should be small. If the conducting 
modes are fully spin effective, i.e., vir ^ wir in Eq. H6.1|l . 
the constant 7 is of the order of the conductivity of the 
junction in units of e^/h. Although smallness of 7 is 
needed for the derivation, the theory is expected to be 
qualitatively applicable even for 7 ~ 1. From this one 
concludes that the intrinsic quantum bound on the mea- 
surement time is given by the inequality, tm ^ fi/eU. 
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get for the lower bound on the measurement time 
AQd (? h 



e 2n\Gi\eU 



(8.7) 



In this case the time of measurement essentially depends 
on the resolution of the charge measuring device. 

In recent experiments involving quantum point con- 
tacts coupled to two-level systems, realized by coupled 
quantum dots, the typical voltage bias was in the mV- 
range and the Rabi frequency of the order of lO^^Hz^^ 
In accordance with the estimate of the measurement 
time, the scheme studied in the present paper can be 
realized in currently studied nanostructures. 



FIG. 3: The conditional entropy for the high voltage case, 
V = 20r, as a function of time for the initial state \^) — 
l/\/2 (!T) + U)), and four different choices of dimensionless 
conductance. The conditional entropy is seen to decrease 
monotonically with a time constant of the order of Tz. The 
entropy curves for gw ~ 0.1, 0.01 and 0.001 all collapse onto 
one curve. 



- - - V/2T=0.1 
V/2T=10 



when the entropy has dropped to one percent, as a func- 
tion of the coupling strength Gw/Gv, for the initial state 
IV') ~ 1/V2(\1) + and high and low bias voltage. It 
is seen that for varying by two orders of magnitude, tm 
roughly scales with Tz as given by Eq. 16.1611 . 



The separation, AQ = \Q^—Qi\, of the two spin orien- 
tations in charge space by the time of measurement, can 
be evaluated as AQ ~ 2\GiU\tm- Using the estimate in 
Eq. iprKjl . 



AQ 

e 



4|Gi 



for any voltage (in the shot noise regime). 

b. A real charge detector, may have a resolution 
lower than AQ given by Eq. (|8.6|l . In this case, the time 
of measurement is estimated from the condition that the 
distance in the charge space between the peaks corre- 
sponding to opposite spin orientations, exceeds the de- 
tector resolution, i.e., AQ{tm) > AQd. From this, we 



IX. SUMMARY AND CONCLUSION 

We have derived the equation of motion, Eq. (|4.5|l . 
for the charge specific density matrix for a quantum sys- 
tem interacting with an environment, the role of which 
is played by a tunnel junction connecting two electron 
reservoirs. In our approach we use the density matrix 
technique with the environment degrees of freedom par- 
tially traced out: We keep track of the charge distribution 
between the reservoirs and derived the master equation 
for the charge specific density matrix of the system, /5„, 
n being the charge variable. The interaction with the 
environment is through the dependence of the tunnelling 
amplitudes on the state of the system. This dependence 
being specified, the derived master equation, Eq. (|4.5|l . 
is general and can be applied to any system. The charge 
distribution is a collective variable of the environment, 
and we have shown that to lowest order in the tunneling 
it can be treated classically. 

In our analysis of the measurement of a spin (two-level 
system), the charge variable plays the role of the pointer 
coordinate in von Neumann's general theory of quantum 
measurementi^ At the conceptual level, the measure- 
ment scheme works as follows: One first disconnects the 
junction, preparing thereby the environment in a certain 
initial charge state, say n = 0. To measure the spin wave 
function \ipo) = a| f) -|- (3\ |) at time t = 0, the spin 
sensitive tunnelling, Eq. H6.1(l . is switched on at f = to- 
gether with the voltage U and kept till the instant t — tm, 
the measurement time discussed in Section rVIIII Then, 
the junction is again disconnected, and the transferred 
charge Q = en, the change in the charge of one of the 
electrodes, is examined. The transferred charge Q is ran- 
dom due to noise, both classical and quantum. However, 
its probability distribution at the instant t = is con- 
centrated in two well-separated peaks located at Q| and 
Ql (as seen in Figures ^(f) and 121(f)). The frequency 
of occurrence of transferred charge Q in the vicinity of 
the corresponding peaks, gives the probabilities jap and 
The measurement is projective: An observation of 
the charge state in a state around Q| (or Q|) ensures 
that the spin after the measurement is in the pure state 
I T) (or I i)). 
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In the shot noise regime, which is the most favorable 
for performing the measurement, the measurement time 
and the voltage U are related as eU tm ^ hf-f where 7 
(see Eq. (IH-SIl ) is the effective coupling constant. For our 
theory to be valid, the coupling constant must be small. 
However, one may expect that the estimate holds even for 
7^1, giving the ultimate quantum bound, eUtm ^ fi, 
for the time duration of the interaction needed in order 
for the tunnel junction to function as a measuring de- 
vice. On the other hand, the duration must be short on 
the scale of the inverse Rabi frequency for the mea- 
surement to be "instantaneous" . For the typical Rabi 
frequencies ~ 10^°iJz,^'^ the measurement becomes 
"instantaneous", that is ^ 1, provided the ampli- 
tude of the voltage pulse U is above the mV level. 

The application of the charge projection method to a 
tunnel junction, a semi-realistic model of a nanodevice, 
allows us to study questions regarding quantum mea- 
surement, and we have considered the emergence of a 
projective spin measurement. We have obtained explicit 
expressions for the characteristic times for decoherence 
and spin state purification, and visualized the emergence 
of a projective spin measurement. The model allowed 
a detailed study of the purification of the charge spe- 
cific spin density matrix. The spin measurement by the 
tunnel junction is quite analogous to the spin measure- 
ment in the Stern-Gerlach experiment, except that the 
interaction mechanism of course is rather different, being 
due to interaction with reservoirs of charged fermions. 
This gives another example of a measuring scheme show- 
ing that there is no need to postulate the classicallity 
of any variable and no "wave function collapse" needs 
to be invoked. The transmission of electrical noise from 
the tunnel junction to the quantum object accomplishes 
the projective measurement. Amplification to the classi- 
cal level and the emergence of a projective measurement, 
have thus been illustrated using a model of a nanodevice. 

The use of charge projectors to study charge kinetics 
originated in the context of counting statisticsiS^ The de- 
rived master equation for the charge specific density ma- 
trix shows that the method of charge projectors leads to a 
useful application of counting statistics. Finally, we note 
that the obtained master equation allows the study of the 
charge specific dynamics of an arbitrary quantum object 
coupled to a tunnel junction. The presented method is 
therefore suitable for attacking problems in a variety of 
fields, such as for example quantum electro-mechanical 
effects in nanostructruresi^ 
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APPENDIX A: DERIVATION OF EQUATION OF 
MOTION 

In this appendix, the equation of motion for the charge 
specific density matrix, /5„, is derived for the case where 
tunnelling events are rare, i.e., to lowest order in the 
tunnelling. It is convenient to use the Heisenberg picture 

p' = e'"^' pe-'"^* (Al) 

where the density matrix evolves only due to the dynam- 
ics of the electrons 

ip'=[H',p'] (A2) 
and the Hamiltonian in the Heisenberg picture is 

H' [t) = Hi + Hr + H'rj,{t) (A3) 
where H!j,{t) = T{t) + T^{t) and 

l,r 

fir[t) = e*-^^*fj^e-*-^s*. (A4) 

The equation of motion in the Heisenberg picture for 
the charge specific density matrix of the system coupled 
to the junction can be written on the form 

Pn{t) = E {TUt)At\t) + f^i{t)B^I^{t)) + h.c. (A5) 

Ir 

where the time dependent operators A and B are 

At\t) = \tt,i {c\crp'{t)rn) , 
Bt\t) = \'^^el {clcip'{t)Vn). (A6) 

A hat indicates operators with respect to the degrees of 
freedom of the system coupled to the junction. 

To obtain, e.g. A, one uses the equation of motion, 
Eq. ||A2|I . to obtain 

Air{t) - lUJlrAirit) = -Tr,j (c^ [H^ (t) , p' (OjT'n) 

(A7) 

where the term containing ljij. = — originates from 
the commutator with the electrode Hamiltonian, Hq = 
Hi + Hr- We must therefore consider the inhomogeneous 
term on the right of Eq. l|A7p . According to the relations 
Eq. I|3.5|l and Eq. (|3.()l) it can be rewritten on the form 

Trei (clcr[{T^ +T}p']Pn') = Trei (clcrTVn-2p'Vn) 
-Tl-el (c\cr.Vn-lp'Vn+ir^) + Tr^i (c\crT^ Vnp' Vn) 
-Tlel (c\crVn^lp'Vn-ir) (A8) 

where the time argument has been suppressed. The first 
two terms are charge "far" off-diagonal components of 
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the density matrix like 'PmP'Pm±2, and these terms can be 
neglected since they are of higher order in the tunnelling 
matrix element and an expression involving only the two 
charge diagonal components of the density matrix results. 
These charge diagonal terms, which have the form 

TTeliclcrT'^Pnp'Vn) = T^/;/ Trej (c^C^ C^, Q/ 7^„p'7^„ ) , 

I'r' 

(A9) 

can be approximated by the expressions 

Tiel {clcrcl,Cl'VnP''Pn) = 5u'5rr'fl{l ~ fr)p'n , 



= Tre/ (VnP'Vn) 



(AlO) 



where the /'s are the single particle energy distribution 
functions for the electrodes. The distribution functions 
are taken independent of the number of charges on the 
left electrode n, since the electrodes are inexhaustible 
particle reservoirs. The diagonality in the electrode quan- 
tum numbers are justified by the fact that when tun- 
nelling is rare, superpositions of different single particle 
states do not occur. 

Approximating similarly the second charge diagonal 
term, Eq. IjASfl . makes Eq. (|A7() acquire the form 

Air - iuJlrAir = - fr)Trl p'n + /r(l " fl)Pn^iTrl- 

(All) 

In the following we shall take as initial condition that 
tunnelling first starts at time t = 0, i.e., the junction is 
disconnected at earlier times. The solution of Eq. (|A11() 
thus becomes 



Air it) = -fl{l-fr) [frl{t)p'Jt) 

+ fr{l-fl) \p'„_,{t)frl{t) 



where the notation 
t 



iujir.(t-t')+i f dt"V{t") 

[X{t)]= / dt'e *' X{t') 



(A12) 



(A13) 



has been used, and the feature that the junction may be 
biased by a time dependent voltage, Vit), is included. 

The evaluation of B in Eq. (|A6p is analogous and we 
obtain (of course in the same approximation and for the 
same initial condition) 



Brl[t) 



fr{l-fl) p'nit)Tlit) 



+ fl{l-fr)\fUt)p'^+,{t)\ (A14) 



where the dagger indicates hermitian conjugation of the 
system operators. 

Collecting the results, we obtain a non-Markovian mas- 
ter equation for the charge specific density matrix: 

lit) = A'{p'^it)}+v'^{p:^^,it)^p:M} 

+ V'_ {p:,_,{t) - Pnit)} (A15) 



where the kernels are 
A' {/5}=E/'(l-/'-)( 



Tlr — Tir 



tIp 



h.c. 



J2fAl-fl) (fir [pfl] - [pfl^fir) + h.c. (A16) 



and 



and 



{P} = E /'(I - [^lrp\^lr + h.c. (A17) 

Ir 



Ir 



h.c. 



(A18) 



The charge specific spin density matrix in the Heisen- 
berg picture evolves with a rate proportional to the elec- 
tron tunnelling rate. The temporal non-locality of the 
kernels in Eq. (|A15|) . however, is independent of tun- 
nelling, and instead depends on the quantum time scale 
determined by temperature and voltage. In the limit 
where tunnelling can be neglected on this time scale, the 
master equation for the spin dynamics becomes Marko- 
vian since p' can be taken outside the bracket operation 
[...]in Eq. (|A16IA18|I . We therefore finally arrive at the 
Markovian master equation, Eq. 14.5|1 . however there 
displayed in the Schrodinger picture. 



APPENDIX B: ONE-COMPONENT MASTER 
EQUATION 



The equations Eqs. (|5.1|l and H6.5I6.8|I . have the struc- 
ture 

Xn = D^x'^ + J^x'^ , a;„(t = 0) = xo (5„,o- (Bl) 
The equation can be solved by Fourier transform. 



Xr,{t) = xo I ^^e(2^-('=°"'^-i)+*'^=='^'"'^)*+*^" (B2) 
In 



givmg 



Xn{t) = Xa 




I„i2t^D,+D,-) (B3) 



where Dx± = Dx ± ^Jx, and /„ is the Bessel function. 

Let Unit) be another variable obeying Eq. IjBljl with 
corresponding coefficients Dy and Jy. Then, the "over- 
lap" of the variables, ^ a:„(t)j/„(t), can be calculated as 

n 

Y.^n{t)yn{t) = xoyoe-2(^-+^«'*/o (f) , (B4) 



where t' = 2t^ {Dx+ + Dy -){Dx - + D y+)- This expres- 
sion is used to derive Eqs. (|6.15|) . and (|6.17|) . 
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1. Charge transfer probability distribution 



as 



The solution to the master equation, Eq. H5.1|l . for the 
probabihty of n-charge transfers, p„, which has the struc- 
ture of Eq. ljBl|) , can be read off Eq. ljB3|) giving 



in the dimensionless variables 



r = 2Dt 



Id 



hi 



D 



(B5) 



(B6) 



With D and / evaluated from Eqs. H5.2|l . and H5.3|l . 

the parameters D and / have the relationship I /2D ~ 
tanh(y/2T), as required by the fluctuation-dissipation 
theorem, and the parameter v has the meaning of the 
dimensionless bias, v = V/2T. 

The charge expectation values, {n^) — ^ n^Pn and 

n 

the moments ((An)*^) where An — n ~ (n), can be 
obtained from the generating functions, and J-u, 

oc oo 

Fu{r) = > ^« = E P»(r)e«("-<"» 

n— — OO n— — oo 

(B7) 



in") 



(B8) 



After summation Eq. ljB7p with j3„ from Eq. (|B5|I , we get 



Fu{t) = exp 



cosh(t; -|- u) 
coshw 



.77 — lirtanhi? 



Fu{r) 
(B9) 



For the first five moments one gets: (n) = 
Ttanhw , ((An)2} = r , ((An)^) = rtanhw , {(An)^) = 
T + 3t^ , ((An)^) = t(1 + IOt) tanhi;. In particular, the 
third moment satisfies the relation 



{{Anf) = {n) 



(BIO) 



in agreement with Ref. |2£ 
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